MIXING TIME BOUNDS FOR ORIENTED KINETICALLY CONSTRAINED SPIN 

MODELS 
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Abstract. We analyze the mixing time of a class of oriented kinetically constrained 
spin models (KCMs) on a d-dimensional lattice of n d sites. A typical example is the 
North-East model, a 0-1 spin system on the two-dimensional integer lattice that evolves 
according to the following rule: whenever a site's southerly and westerly nearest neigh- 
bours have spin 0, with rate one it resets its own spin by tossing a p-coin, at all other 
times its spin remains frozen. Such models are very popular in statistical physics be- 
cause, in spite of their simplicity, they display some of the key features of the dynamics 
of real glasses. We prove that the mixing time is O(nlogn) whenever the relaxation 
time is O(l). Our study was motivated by the "shape" conjecture put forward by G. 
Kordzakhia and S.P. Lalley. 
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1. Introduction 

Kinetically constrained spin models (KCMs) are interacting 0-1 particle systems, on 
general graphs, which evolve with a simple Glauber dynamics described as follows. At 
every site x the system tries to update the occupancy variable (or spin) at x to the 
value 1 or with probability p and q respectively. However the update at x is accepted 
only if the current local configuration satisfies a certain contraint, hence the models 
are "kinetically constrained". It is always assumed that the constraint at site x does 
not to depend on the spin at x and therefore the product Bernoulli(p) measure tt is 
the reversible measure. Constraints may require, for example, that a certain number 
of the neighbouring spins are in state 0, or more restrictively that certain preassigned 
neighbouring spins are in state (e.g. the children of x when the underlying graph is 
a rooted tree). 

The main interest in the physical literature for KCMs (see e.g. [20] for a review) 
stems from the fact that they display many key dynamical features of real glassy ma- 
terials: ergodicity breaking transition at some critical value q c , huge relaxation time 
for q close to q c , dynamic heterogeneity (non-trivial spatio-temporal fluctuations of 
the local relaxation to equilibrium) and aging, just to mention a few. Mathematically, 
despite their simple definition, KCMs pose very challenging and interesting problems 
because of the hardness of the constraint, with ramifications towards bootstrap percola- 
tion problems [23], combinatorics [7, 24], coalescence processes [10, 11] and random 
walks on upper triangular matrices [18]. Some of the mathematical tools developed 
for the analysis of the relaxation process of KCMs [4] proved to be quite powerful also 
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in other contexts such as card shuffling problems [2] and random evolution of surfaces 
[5]. 

In this paper we focus on oriented KCMs on a ci-dimensional lattice, d > 2, with n d 
sites, in particular on their mixing time. A prototypical model belonging to the above 
class of KCMs is the North-East model in two dimensions (see e.g. [12] and [4]) for 
which the constraint at any given site x requires the south and west neighbours of x to 
be empty in order for a flip at x to occur. In order to avoid trivial irreducibility issues 
the south-westerly most spin is unconstrained and sites outside the upper quadrant are 
treated as fixed zeros. 

With p c the percolation threshold for oriented percolation in two dimensions (see 
e.g. [9]), it was proved in [4] that for all p < p c the relaxation time of the North-East 
process is O(l) while it becomes life™) 1 , for some c > 0, when p > p c . Atp = p c the 
relaxation time should have a poly(n) growth. Consider now the North-East model in 
the first quadrant of Z 2 . In [12] it was conjectured that, for p < p c and starting from all 
l's, the influence region R t , defined as the union of all unit squares around those sites 
which have flipped at least once by time t, has a definite limiting shape S C M? in the 
sense that a.s. ^ — > S as t — > oo. Since the North-East process is neither monotone or 
additive (see [14]), the usual tools to prove a shape theorem do not apply in this case. 

The above conjecture implies that, for p < p c , the influence coming from the uncon- 
strained spin at the South-West corner propagates at a definite linear rate as it does 
in the East model [3], the one dimensional analog of the model (for background see 
[1, 4, 22]). In particular the mixing time of the model should grow linearly in n (the 
linear size of the system) . However in dimension d > 2 the analysis of the propagation 
of influence is quite delicate because of the many paths along which it can occur (see 
[24] for combinatorial results in this direction) . 

In this paper we prove that the mixing time is 0(n log n) as long as the spectral gap 
of the process is 0(1). Our technique bares some similarities to those employed in [6] 
to analyse the Glauber dynamics of biased plane partitions. 

2. Models and Results 

2.1. Setting and notation. We consider a class of 0-1 interacting particle systems on 
finite subsets A of the integer lattice Z d , reversible with respect to the product measure 
7r := n^eA 7T ' X ' wri ere ir x is the Bernoulli(p) measure. 

The (f-dimensional cube of linear size n (which contain n d points) will be denoted 

by 

A n : = ([l,n] x...x[l,n])nZ d . 

The standard basis vectors in Z d are denoted ei = (1, . . . , 0), ei = (0, 1, ... , 0), . . ., 
ed = (0, ... , 1). For x G A n we write Xj for the component of x in the direction e,. 

The set of probability measures on the finite state space Q n = {0, l} An is denoted by 
V(Q n ). Elements of f2 n will be denoted by the small greek letters a,r/, . . . and a x will 
denote the spin at the vertex x. 

We denote the i-th level hyperplane in A n by Hi = {x € A n : Y^=i x j = ( see 
Fig. 1). The set of sites on and below this hyperplane will be written Ui = {x <E A n : 

1 We recall that / = Q(g) if |/| > eg for some c > 0. 
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Figure 1. Left: The the lattice A n in two dimensions, d = 2. The 
sets Hi and Wj-i are shown in dark grey and light grey respectively. 
Right: The constraining neighbourhoods, fC x and K.*, of a vertex x in 
two dimensions, shown in dark grey. 



Ylj=i x j ^ ^ or eacri c G f2 n we write o-W for the restriction of <r to Ui, and for 
the restriction of a to Hi. Similarly, for any probability measure v G V(£l n ), we write 
i/W for the marginal of on : = {0, 1}"* c O n . 
For any vertex x e A n we also let (see Fig. 1 Right) 

/C* = [y G Z d : y = x - ^a.e,, a, G {1,0}} \ {x} 

i=l 

K-x = \ y G % d '■ 3 i = 1, . . . , d such that y = x — eA . 



Definition 2.1 (constraints). Let x* = (I, ... ,1) be the south-west corner of A n . Consider 
a collection {C x } xe \ n of constraining neighborhoods such that 

C*. = 0, and <D^C x cK*n A n , x ^ x*. 

Let 

We then say that the constraint at site x is satisfied by the configuration a if c x (a) = 
1. In words, the constraint at x is satisfied if all the spin in C x are 0. Note that x* is 
unconstrained. 
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2.2. Oriented KCMs and main result. We give a general definition of the process to 
include a large class of directed KCMs, such as the North-East model and higher dimen- 
sional analogues. For constraining neighborhoods {C x } xe \ n we define the associated 
directed KCM by the following graphical construction. To each x G A n we associate 
a mean one Poisson process and, independently, a family of independent Bernoulli(p) 
random variables {s x ^ ■ k G N}. The occurrences of the Poisson process associated 
to x will be denoted by {t Xy k '■ k G N}. We assume independence as x varies in A n . 
With probability one all the occurrences {t Xi k}ken,x&\ n are different. This defines the 
probability space and the corresponding probability measure is denoted by P. 

Given r\ G f2 n we construct a continuous time Markov chain {rj(s)} s j> o on the prob- 
ability space above, starting from rj at t = 0, according to the following rules. At each 
time t x>n the site x queries the state of its own constraint c x (see (2.1)). If and only if 
the constraint is satisfied (c x = 1) then t x>n is called a legal ring and the configuration 
resets its value at site x to the value of the corresponding Bernoulli variable s Xjn . 

The above construction gives rise to a irreducible, continuous time Markov chain, 
reversible w.r.t. tt, with generator 



where n x denotes the conditional mean 7r(/ | {<T y } y ^ x ). Irreducibility follows because 
we can invade A n with 0's starting from the unconstrained corner x*. 

Using a standard percolation argument [15, 8] together with the fact that the con- 
straints {c x } x6 a„ are uniformly bounded and of finite range, it is not difficult to see 
that the graphical construction can be extended without problems also to the infinite 
volume case. 

For an initial distribution u at t = the law and expectation of the process will be 
denoted by ¥ u and E„ respectively. In the sequel, we will write v t for the distribution 
of the chain at time t, 



If v is concentrated on a single configuration ?] we will write P^-)- Note that, for each 
i = 1 , . . . , n, the same graphical construction can be used to define the process on 
whose law is denoted by (■). 

It follows from the graphical construction, and the fact that the constraints are ori- 
ented, that given an i G {1 . . . , n} the evolution in Ui is not influenced by the evolution 
above Ui. In particular, for any i] € Vt n and any event A in the a-algebra generated by 

{nx(s)} s ^ t, x <=Ui, 



In fact the same holds for any subset U C A„ with monotone surface, i.e. x — G U 
whenever x G U and x — G A n , for each i G {1, . . . , d}. 

We finish this section with definitions of the spectral gap and mixing time of the 
process. 




(2.2) 



v t [rj\ :=¥ u {i 1 {t)=i 1 ). 



(2.3) 
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Definition 2.2 (spectral gap). The spectral gap, A^ n ^ of the infinitesimal generator (2.2) 
is the smallest positive eigenvalue of —C n , and is given by the variational principle 

AW:= inf (2.4) 
/:Q„^R Var n (/) 

fyt const 

where V n (f) = —tt (f£ n f) is the Dirichlet form of the process. 

Definition 2.3 (Mixing time). The mixing time is defined in the usual way as 

T £L = > : SU P IK " tIItv < V4} , (2.5) 
uev{n n ) 

where \\ ■ ||tv denotes the total variation distance. 

Remark 2.4. It follows from [4, Theorem 4.1] that there exists < po < 1 such that 
inf n A^ n ) > whenever p < pq. For the North-East model po coincides with the oriented 
percolation threshold. Moreover, for any p <G (0,1), one has l/A^ ^ ^ cn d /X^ 
(see e.g. [13]j. 

Our main result reads as follows. 

Theorem 2.5. Assume inf„ A^ n ) > 0, then 

Cn^T^l^C'nlogn (2.6) 
for some C, C > independent of n. 

Remark 2.6. 7n the case of maximal constraints, i.e. C x = K* x fl A n , inserting the indicator 
of the configuration identically equal to 1 as a test function in the logarithmic Sobolev 
inequality [21] shows that the logarithmic Sobolev constant c/ n ) of the generator (2.2) is 
0(l/n d ). Therefore, if 

T 2 (n) := m£{t > : supVar ( — ) < 1/4}, 
then for a universal constants c > and a constant d > depending only on p, 

where the first inequality follows from [21, Corollary 2.2.7]. This observation shows that 
oriented models in Z d can be quite different from oriented models on rooted k-regular 
trees. In the latter case it was recently proved [16] that both the mixing time T ^v and 
grow linearly in the depth n of the tree whenever the relaxation time 1/A^ is O(l). 

3. Proof of theorem 2.5 

Proof of the lower bound. The lower bound on the mixing time is quite straightforward 
using finite speed of propagation. Define r* as the first time the spin at site x* := 
(n, . . . , n) is 0, and denote the configuration of all l's by 1. Using results in [17, 19] 
there exists c > such that 

T !nl > cEi [t*] > (r* > n/2) . (3.1) 
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Clearly the event r* < n/2 requires the existence of a path 7 = {x* = , , . . . , x^> = 
x*} and times to < t\ < ■ ■ ■ < t£ < n/2 such that, for any i ^ £, xW € C z (i+i) and ij is 
a legal ring for a;W. Standard Poisson large deviations show that the probability of the 
above event is o(l) as n — > 00. □ 

Proof of the upper bound. The proof of the upper bound is based on an iterative scheme. 
For i < n, we find some time, Aj = 0(log(z)), such that if the initial measure v has 
marginal equal to vr^" 1 ) on VL^ ^ (i.e. below the hyperplane Hi), then after time 
Aj the marginal i/W is very close to the equilibrium marginal tt^\ This is the content of 
Lemma 3.1. Then, starting from an arbitrary initial measure, we can iterate the above 
result using the triangle inequality for the variation distance and propagate the error. 

Lemma 3.1 (Mixing time for a single diagonal.). There exists c = c(q, d) > such that, 
for any initial measure v with marginal on O^ 1 ) equal to / k^ % ~ 1 \ 



TV 



$C e for all t c log(z/e) . 



Before proving Lemma 3.1 let us recall a useful characterisation of the total variation 
distance (see for example [13]), 



W- HItv = o su p K/) _7r (/)l> 

2 ll/lloc < 1 



where ||/||oo denotes the sup-norm. 

Proof of Lemma 3.1. Fix i ^ n and a function / depending only on the spin configu- 
ration in Ui and such that ||/||oo ^ L Without loss of generality assume tt(/) = 0. 
Given an initial measure ^ with marginal vr^" 1 ) on f^" 1 ) it follows from (2.3) that 

f t = 7T^ _1 ) for all t ^ 0. Moreover, conditioned on the history {cr( l_1 )(s)} s ^ t 
in the spins {cr x (i)} xe ^ i on the hyperplane % evolve independently from each 
other. Each one goes to equilibrium with rate one during the intervals of time in which 
its constraint c x is satisfied and stays fixed otherwise. In particular, conditioned on 
having had a legal ring at each x e H% before time t, the distribution of a-^ i (t) is it. 
These simple observations gives rise to an upper-bound on the expectation of / at time 
t as follows. 

Let rW be the first time there has been at least one legal ring on each x € Hi. 
Following the argument above we have, 



+ 



K (E 
E„ (E 
E„(E 



/(>(*)) l {rW<t} 



s < t 



E, E 



{rW > t} I {^(s)} 



+ 



L {rW>t} 



< E, 



In the second inequality above we used the strong Markov property together with the 
observation that the distribution of o?^(£) conditioned on {rW < t} and {^^^(s)} s < 
is 7r and the fact that 7r(/) = 0. 

To bound the final term above we denote the number of rings of the Poisson clock at 
site x during the set B c [0, t] by N X (B) = \{t Xt k}k > o n B\, and we define the set of 
legal times at site x before t as Q(x,t) := {s < t : c x (a(s)) = 1}. By construction, for 
x E Hi the set G(x, t) depends only on {u y (s)} s <t , y G C x . Using the notation above, 

K^ytHa^Hs)}^,) < ^p(iV a .(e(x, i ))=0|{^- 1 )( S )}^ i ). 

By construction N x (Q(x, t)) is a Poisson random variable of mean \<3(x, t)\. Thus 

\S(x,t)\ 



Y J V(N x (g(x,t))=0\{a^ 1 Hs)} 

xeHt 



s < t 



x£Hi 



In conclusion, 



E„ 



<; z d 1 max E w 
yeA n 



\S(y,t)\ 



We estimate [e I^Cs/>*)n using a Feynman-Kac approach (a similar method has been 
used to bound the persistence function in [4]). The total time that site y satisfies its 
constraints is 



\G(y,t) 



[ c y (a(s)) 
Jo 



ds. 



On L 2 (ir) we define the self adjoint operator H := C — V with V(a) := c y (a). The 



as 



<1> 



e tH l) 



Feynman-Kac formula allows us to rewrite the expectation E n 

(where (•, -) n denotes the inner-product in L 2 (ir)). Thus, if j3 y is the supremum of the 
spectrum of H we have 



E, 



\G(y,t)\ 



(3.2) 



In order to complete the proof of the Lemma it remains to show that f} y < 0. 

We decompose each L 2 (7r)-norm one function <f> in the domain of H as <f> = al + g, 
for some mean zero function g. So (1, g) n = and a 2 + (g, g) n = 1. Then, 

(</>, H<P) = (g, Cg) n - a 2 (1, VI) n - 2a (1, Vg) w - (g, Vg) n 



(3.3) 



We proceed by bounding {<p,H(p) n . Using (3.3), 

(<P, H<P) < - (g, g)^ - q d (a 2 + 2vr (g \ V = 1) a + vr (g 2 | V = l) 
5X (n \ 



where 5 = 1 — a 2 . Above, in the second inequality we used the fact that the term in 
the large brackets (with common factor q d ) is a quadratic over a which is non-negative 
when q = and has non-positive determinant by Jensen's inequality. Again from 
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(3.3), dropping the last term on the right hand side and applying the Cauchy-Schwarze 
inequality we get 



H4>) < - A< n > (g, ~ "V + 2|a| ((s, g),q d (l ~ q d )) 
< (6 - l)q d + 2 f <5^(1 - £ 



1/2 



1/2 



Thus 



Py *S - c (g) 



mm 



-<5A (n) ,(<5-l)^ + 2(V(i-g d ) 



1/2 



Finally we observe that co(q) > because 



(5-1)/ + 2 
In conclusion 



1/2 



if J< 



2(1 - ^l-g Q 



4g a 



E, 



/(' W (t) 



for i ^ clog(z/e) for some c = c(d,q). 



□ 



We now prove Theorem 2.5 by iterating the previous Lemma and propagating the 
error. We use the following iterative scheme 

t i = t i _ 1 + A i , t = 0, (3.4) 

where Aj(e) = Mog(i 3 /e) with c as in the Lemma. We now show by induction that 
starting from an arbitrary initial distribution v 



TV 



* 1 

^ ii := e ^ f° r a ^ t ^ *i 



(3.5) 



.7 = 1 



The case i = 1 (3.5) is an immediate consequence of the fact that the corner x* = 
(1, . . . , 1) goes to equilibrium with rate one. Assume (3.5) holds for all i ^ k — 1. It is 
clear that min CT ft (a) > for any t > 0, so that we may define an auxiliary probability 
measure fj, with marginal ^ k ~ l "> = tt^" 1 ^ by 



7T 



(fc-1) 



(7 



(fc-1) 



"t*-l 



0" 0" 



(fc-1) 



Fix a function / depending only on the configuration in U\~ with 
without loss of generality, assume ir(f) = 0. Then 



E„ [/ (<r(t fc ))] 



E, 



[/(*(**))] 



E^ [/(a(A fe ))]-E^j>(A fc ))] 



+ 



A- 2 



< 1. Again, 



(3.6) 



where we used the Markov property in the first line and Lemma 3.1 together with the 
triangle inequality in the second line. The two measures // and ^ fc _ 1 have the same 
marginal on Hk, so we may reduce the remaining term to something that can be dealt 
with using the inductive hypothesis as follows. For any 77 € fj( fc_1 ) let 



F(n) :=E, 



f(a(A k )) \a^= V 
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Clearly ||.F||oo ^ 1, so by the definition of fi and the inductive hypothesis (3.5), 

[/ (a(A k ))] - E„ [/ (a(A k ))] | = ^ (F) - tt (F)\ ^ 4_i • 
In conclusion the r.h.s. of (3.6) is bounded from above by 

K [/(^(tfc))] 

Since Ui n -\ 5 A n , it follows that 



2n-l 



2 ||ft - 7r|| TV < e 2n _i = e < 2e for all i ^ i 2ri _i 



and 



hn-x = y~) Aj - / log(x 3 /e)(ix < c'nlog(n) 
i=i C - 70 



for some c' > 0. 



□ 
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